Directed diffusion of reconstituting dimers 
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Q. 

We discuss dynamical aspects of an asymmetric version of assisted diffusion of hard core particles 
on a ring studied by G. I. Menon et al. in J. Stat Phys. 86, f237 (1997). The asymmetry brings in 
phenomena like kinematic waves and effects of the Kardar-Parisi-Zhang nonlinearity, which combine 
with the feature of strongly broken ergodicity, a characteristic of the model. A central role is played 
by a single nonlocal invariant, the irreducible string, whose interplay with the driven motion of 
. reconstituting dimers, arising from the assisted hopping, determines the asymptotic dynamics and 

scaling regimes. These are investigated both analytically and numerically through sector-dependent 
mappings to the asymmetric simple exclusion process. 

PACS numbers: 05.10.Gg, 02.50.-r, 05.40.-a, 02.50.Ey 
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The issue of universality classes in noncquilibrium statistical systems is often linked to the existence of conservation 
laws Normally, the number of conservation laws is finite, leading to the occurrence of dynamically disjoint 
sectors whose number grows as a power of the system size. However, certain dynamical processes involving composite 
objects exhibit strongly broken ergodicity, with the number of disjoint sectors growing exponentially with system 
size, as a result of having an extensive number of conservation laws. Examples studied earlier include deposition 
and evaporation (•••<-+ooo) 2,3.4.5], and diffusion (• • o •(->• o • •) Q of bunches of particles. If these moves are 
interpreted as involving trimers or dimers (in general fc-mers), then the fc-mers in question do not keep their identity, 
and can reconstitute in time. The dynamical moves in these models do not connect configurations in different sectors, 
implying that the steady state is not unique and depends strongly on the initial condition. Moreover, the form of 
the long-time decay of time-dependent correlations in the steady state varies strongly from one sector to another. In 
one dimension, both the partitioning of phase space into many sectors and the accompanying dynamical diversity 
could be understood in terms of a nonlocal construct known as the Irreducible String (IS) which is an invariant of 
^sO [ the motion. The IS provided a convenient label for each sector Moreover, the position of the elements of the 
IS were the relevant 'slow' variables in the problem, and thus governed the long time dynamics in different sectors. 
The dynamical diversity found in different sectors (ranging from different power law decays to stretched exponentials) 
could be accounted for in terms of the differences in the IS from sector to sector [j| . In turn, the IS can be used to 
construct an extensive number of conservation laws, although it turns out that these involve nonlocal combinations 
of the site occupancies Q. 

All such studies known to us involve a symmetric movement of the IS. Here we study the dynamical consequences 
of an asymmetric (directed) motion of the IS. This brings in new phenomena associated with driven diffusive systems, 
O ■ such as kinematic waves and effects of the Kardar Parisi Zhang (KPZ)/Burgcrs nonlinearity, present in the 1-d 
asymmetric simple exclusion process (ASEP) 0§|. We investigate these effects by studying the directed diffusion of 
reconstituting dimers (DDRD). The model resembles the diffusing reconstituting dimer model studied in Ref. [6|, the 
only difference being that we allow for only forward motion of the dimers (• • o — > o • •), in contrast to the two-way 
motion studied in that paper. The construction of the IS is the same for symmetric/asymmetric motion, but the 
dynamics of dimers (and of the IS) is quite different in the directed case. 

It turns out that there is a correspondence between the DDRD and the well-studied ASEP, though of a generalised 
sort, leading to extra features including additional wheeling motion of sites, oscillations of correlation functions, and 
other new kinetic effects. These typically result from the combination of collective driving and the invariant, but 
moving, IS. Interestingly, there is a very useful correspondence between the DDRD and a three-species exclusion pro- 
cess. Further, the problem of forward-moving non-reconstituting (hard) dimers which keep their identity is recovered 
in the DDRD, in a particular sector. 

As in the previous studies of dynamic diversity, Monte Carlo (MC) simulation is a key ingredient in trying to 
understand the dynamics. In particular we study time-dependent correlation functions in different sectors and use 
analytic reductions related to the specific form of the IS and the motion of its elements to correlate the diverse 
behaviour seen in the simulations with properties of the ASEP. Remarkably close correspondences in nonuniversal as 
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well as scaling properties are seen. 

The development begins (Section 2) with an introduction to the model and then moves on in the same Section 
to the following sequence of emerging topics: equivalence to the three-species process and correspondences with the 
ASEP, whose properties are summarized. In Section 3 we study the kinetics within a particularly simple sector 
in which the problem is tantamount to that of non-reconstituting dimers. New effects such as wheeling of sites in 
the equivalent ASEP, resulting from the combination of driving and the composition of dimers, are shown to have 
important consequences for autocorrelation functions. These show early oscillations and later a decay, whose form 
(exponential or power law), is decided by a critical condition related to wheeling and kinematic wave velocities. The 
long time behaviour shows scaling, and universality in the sense of data collapse to ASEP scaling functions, but with 
sector-specific parameters. The theoretical and MC results are extended to more general sectors in Section 4, via 
investigations of sublattice currents, sublattice current-density relations, and consequent kinematic wave velocities. 
Section 5 is mainly concerned with spatial correlations in a particular sector, involving MC results and their exact 
analysis via one of the equivalent models. Section 6 is a concluding discussion. 



The DDRD model consists of a ring of L sites, each of which may be singly occupied (occupation variable rij = 1) or 
empty (rii — 0). Any of the 2 L possible configurations is then an L-bit binary string. The system evolves stochastically 
through the move 110 — * Oil, which represents the directed diffusion (rightward hopping), of reconstituting dimers 
(DDRD). This is equivalent to stochastic hopping of holes two steps to the left (i.e. staying on the same sublattice) 
provided the intervening site is occupied. Theprocess is the totally asymmetric alternative to the fully symmetric 
reconstituting dimer diffusion process (DRD) |6| which is known to be strongly nonergodic due to the existence of a 
conserved string of variables, the IS. All states linked by the dynamic process have the same IS, so the phase space 
divides into many sectors, ~ X L in number, where A is the Golden number (vo+ l)/2 . All these properties are shared 
by the asymmetric generalisations. The IS for a given sector can be obtained from any state in the sector by deleting 
from its L-bit string any pair of adjacent l's, and repeating the procedure until no more deletions are possible (the 
result is independent of the order of deletion). For instance, the configuration C = 11101001111010 leads to the IS 



An equivalent representation of configurations ,of the IS, and of the process, uses characters A, B, C related to 
the binary variables by A = 11, B = 10, C = [6]. With the periodic boundary conditions used here any DDRD 
configuration can be uniquely decomposed into a configuration of j4's, -B's, and C's. For example, the configuration 
C of the last paragraph can be written as ABBCAACB. Moreover, the totally asymmetric dimer hopping move 
110 — > 011 corresponds to either AB — > BA, or AC — > CA; the former move involves dimer reconstitution. The IS 
construction now corresponds to deletion of all A's, so the IS is a string of only B's and C's, e.g. for the configuration 
C defined above the IS is BBCCB. The absence of any exchange of B's and C's in the DDRD process verifies the 
conservation of the IS. 

In terms of the characters A,B,C, the DDRD process has an obvious and important analogy with the (totally) 
asymmetric exclusion process (ASEP), in which mutually excluding particles hop to nearest neighbour vacancies on 
the right. In this correspondence the dimers A play the role of the ASEP particle, and the B's, and C's of the IS 
correspond to ASEP vacancies. In order to exploit the correspondence here and later, we use Na, N b , Nc to denote 
the numbers of dimers, 10 pairs, and single 0's, respectively. Then the lengths of the DDRD lattice, the IS, and the 
ASEP lattice are L = 2N A + 2N B + N C , £ = 2N B + N c , and L x = N A + N B + N c , respectively. 

The particle density in the ASEP is x = Na/Lx, and the fraction of zeros of C type is y — Nc/(Nb + Nc)- Then, 
in the case of a periodic IS, if the periodic unit [...] of the IS contains n B B's and nc C's, we have y = nc/{n B +nc)- 
It is straightforward to show from these definitions that L/Lx = 2 — y(l — x). This ratio provides a metric factor 
particularly important for converting known currents, tagged hole velocities, and particle densities for the ASEP into 
corresponding DDRD quantities. An important point is that the ASEP image of a fixed DDRD site is not fixed, but 
moves forward and wheels around the periodic ring. The causes and consequences of wheeling are discussed in detail 
in Section 3. 

We conclude the section by recalling some facts about the ASEP that we will need later. On a ring, the ASEP 
has a product-measure uniform steady state (SS), with current Jx = x(l — x) at ASEP particle density x. Density 
fluctuations move as a kinematic wave through the system [8.9], with velocity U = d x Jx = (1 — 2x). Also, in the long 
time scaling regime, the ASEP density-density correlation function Cx(r,t) = (n(r, t) n(0, 0)) — (n) 2 is of the form 



II. DDRD MODEL AND CORRESPONDENCE TO ASEP 
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It is known [Eq. (4.8) of [T^]] that in the limit of large s, g"(s) ~ exp ( — [i |s 3 |) with /x ~ —0.295 , so the correlation 
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function is exponentially decaying in time at large time t. The exception is when the co-moving condition r = Ut 
applies, a special case being the autocorrelation function (r = 0) when x = 1/2 (making U — 0). In these cases the 
correlation function shows i -2 / 3 power law decay in time. 

In view of the DDRD-ASEP correspondence, one might expect similar behaviour for asymmetrically diffusing dimers, 
unlike symmetrically diffusing dimers of Ref. [||, where power law decay is always seen. Detailed considerations of 
the correspondence show that the results are broadly in agreement with these expectations, but modified by a variety 
of interesting extra features discussed below. 

III. NULL SECTOR 

We denote a sector whose IS contains no B elements as a null sector. In terms of the notation introduced earlier, 
this refers to [0] c with C = Nc and y = 1. In any such sector, the DDRD with an even number of particles is 
isomorphic to directed diffusion of nonreconstituting or 'hard' dimers. The reason is that every cluster of l's must 
contain an even number of particles (otherwise the IS would not be all 0's), and hence a cluster can be labeled as 
DEDEDEDE, with a I? at the start and E at the end. Each successive DE pair can be thought of as a dimer, 
which retains its identity forever: the elementary move 110 — > Oil is DEO — ► ODE. The problem of nonreconstituting 
hard dimers (and more generally fc-mers) is of interest in several contexts, and some aspects have been studied earlier 
[lolflll ]. This isomorphism allows us to obtain the correlation function for the reconstituting dimer problem, in the 
special case of the null sector, from the derivation given in Section 5 for the non-reconstituting case. 

We now summarize the results for dynamic properties of interest in null sectors. A more detailed account will 
appear in [jjjj . By following the motion of a hole in the DDRD and its image in the ASEP, we find that the ratio of 
the currents J DDRD (for DDRD), and Jx is equal to the metric factor Lx/L = 1/(1 + x). Evidently, we then have 

J DDRD =X(1-X)/{1+X). 

Wheeling effect. An important effect shows up in a more detailed consideration of the mapping between sites in 
the ASEP and in the DDRD. A fixed site in the DDRD problem corresponds to a moving site in the ASEP - it 
wheels around the ring at a finite mean velocity. This happens because if we follow the ASEP image of a single jump 
110 — > 011 in the DDRD, we find that the image of the central site of the triplet advances by one unit, while the images 
of the other two sites remain unchanged. In time t, the displacement of a mapped site is thus Ar(t) — Wt + <fi(t) 
where W is the wheeling velocity and <p(t) is a zero-mean variable which represents the effect of stochasticity in the 
motion. It can be shown [l3T ] that in the null sector, the wheeling velocity is given by 

W = x(l — x)/(l + x) . (2) 

The wheeling motion of ASEP sites has been verified by direct observation of mapped site motion in MC simulations, 
which shows a small jitter corresponding to cj>(t) around the predicted average wheeling velocity. 

Autocorrelation functions. The considerations of the previous paragraph have important consequences for the 
DDRD: correlation functions are like those of the ASEP only if wheeling of sites (r — > r + Wt) is allowed for. This 
suggests that in the long time scaling regime, if the effects of jitter arising from the stochastic part <fi(t) can be ignored, 
the density-density correlation function C (r, f) of the DDRD in the null sector will be of the form 

C(r, t) = r 2 l z F( [r + (W-U)t] r 2 / z ) . (3) 

For the autocorrelation function (r — 0) the argument of F is proportional to the difference of the wheeling velocity 
and the kinematic wave speed U of the ASEP, W — U = x(l — x)/(l + x) — (1 — 2x) . Then, if F is similar to the 
scaling function g" of the ASEP, the autocorrelation function C(r = 0,t) = C(t) of the null sector will decay as 
i -2 / 3 exp(— nt) at large times, where n is a constant proportional to W — U. On the other hand, in the special case 
where W — U vanishes, late-time slower-than-exponential decay should be expected, that is C(t) cx t~ 2 ^ 3 Y(t) where 
Y falls slowly, e.g. approaching a constant. Notice that W — U vanishes at the 'compensating' value x c of the ASEP 
particle concentration given by 

x c = V2-l, (4) 

which would correspond to an IS [0] £ of length C — L (y/2 — 1) . 

Numerical simulations . To support these expectations, we have conducted extensive simulations of forward- 
diffusing, reconstituting dimers in a variety of situations. Since autocorrelation functions of interest refer to sponta- 
neous fluctuations in steady state (SS), we avoided the relaxation of the system by generating directly SS configu- 
rations, thus saving considerable CPU time. In the equivalent ASEP, all configurations in SS are equally weighted 
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under periodic boundary conditions and can be translated to DDRD configurations. This enables us to average our 
measurements over a large number of samples to reduce scatter in the data. 

In null sectors with IS of length C, such initial conditions are obtained by random deposition of (L — C)/2 monomers 
on a ring of [L + C)/2 sites. Subsequently, each monomer is duplicated by adding another one over an extra adjacent 
location specially added for that purpose, i.e. 1 — > A = 11 , which leave us with a ring of L sites and (L—C)/2 randomly 
distributed dimers. Notice that this generating procedure is not equivalent to a random sequential adsorption of dimers 
in the original system. Had the latter method been applied it would necessarily introduce correlations between dimers 
because of shielding effects among them [l4lfl5| . 

We then evolved our DDRD system using the stochastic microscopic rules referred to above. After a sequence of 
L update attempts at random locations, the time scale is increased by one unit, i.e. t — > t + 1 , irrespective of these 
attempts being successful. Also, measurements at shorter (longer) time intervals can also be allowed by updating 
the system in M < L (M > L) microsteps, and increasing t as t + M/L. Typically, we considered rings of 10 5 
sites and averaged our measurements over 2 x 10 4 histories starting from the independent configurations constructed 
above. In Fig.Q]we show the resulting autocorrelation functions for various particle densities p in null sectors of length 
C = L (1 — p) , or equivalently, at different concentrations x = p/(2 — p) of the associated ASEP system. The full line 
corresponds to the compensating condition x c = \[2 — 1 , thus making U — W = 0, and is the case where slower than 
exponential decay is seen. Specifically, it is consistent with a large time behavior C{t) oc t~ 2 / 3 , in turn suggesting 
that the slow decaying function Y(t) referred to after Eq. (3) approaches a constant at large t. This is in marked 
contrast to all other IS fractions studied which always give rise to exponential decays, either below or above x c . This 
constitutes a strong confirmation of the theory. 

In preparation for the analysis of SS currents in more general situations (see Section 4) , we also estimated numerically 
these quantities by measuring the three-point correlations involved in the assisted hopping, namely 

2 2 

Jl = T )^ (n 2 i-l 7l2i«2i+l) , Jl = £ (n 2t Il2,+l«2j+2) , (5) 

i i 

on each of sublattices 1 and 2, where n = 1 — n denotes a vacancy. As will be discussed below, the compensating 
condition also corresponds to the vanishing of both current derivatives with respect to particle densities. The inset of 
Fig. Q] displays the currents of sector [0] (both equivalent in this case), as a function of the particle density {p\ = p 2 ), 
closely following the J DDRD current discussed above in this Section. As expected, its extremum occurs very near to 
p c = 2 — \/2 namely, close to x — x c . 

DDRD and ASEP analogies. The above results would indicate that the long time autocorrelations of the DDRD 
and ASEP systems are similar provided that the kinematic wave in each system is stationary. Despite the existence 
of jitter between site locations in the two systems, Fig. [5] provides a strong indication that this might be the case 
provided large times are considered (t > 20) . On the other hand, the short-time differences between both systems (see 
inset of Fig. [2] ) can be qualitatively understood by combining wheeling with the equal-time spatial structure, which, 
as will be discussed in Section 5, has a finite SS correlation length. Moreover, analogues with the ASEP would suggest 
attempting a collapse of all the null sector DDRD autocorrelation data at x ^ x c . Fig. [3] shows the data collapse for 
the DDRD autocorrelation C(t) at representative x's in the null sector, as a plot of \n[bt 2 ' 3 C(t) ] versus a t. There is 
a convincing collapse to a line corresponding to ln<7"(s) versus s 3 , where g" is the ASEP scaling function [Eq . (fTJ)] in 
which s 3 oc t . As for a and b, they are fitti ng p arameters which depend on x, and their numerical values [13l ] closely 
follow those arising from the analysis of jl2lil6l ]. 

The argument s of g"(s) and the overall prefactor required to give the autocorrelation function for the ASEP are 
given on dimensional grounds as functions of x in [l2j and [l6[ (in particular, see Eqs. (1.3) to (1.10) in [16] ), and 
involve U = (1 — 2x) and Jx = x(l — x) together with t in forms consistent with the scaling variables given in Section 
2. There are also numerical factors, which differ between the two papers. Our data collapse for the DDRD yields 
values for a which for most x are within 1 per cent of the corresponding theoretical ASEP values [16fl ; likewise for 6, 
except for a persistent difference of a factor of 2 (possibly owing to the definition of g" in [16]). Apart from this, the 
DDRD autocorrelations correspond in every detail to the ASEP ones except at early times. 

The data in Fig.[3]do not lie in the asymptotic regime where In g"{s) is linear in s 3 . This exponential decay is achieved 
only for s 3 > 200, as displayed in upper inset; the corresponding times are very large and the autocorrelations too 
small as to render them numerically accessible by direct simulation. Size dependences are then noticeable unless 
L 3 / 2 is large compared to such times, so very large system sizes are required in the MC. Further, the implied large 
equilibration times mean that for large L it would not be feasible to prepare initial states by relaxation. In our 
simulations, initial states are obtained via the ASEP by exploiting its product measure property, as described above. 
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IV. PERIODIC IRREDUCIBLE STRINGS 

In non-null sectors, the value of Nb is a measure of reconstitution that occurs in that sector. We have examined 
some sectors with periodic IS's and find pronounced oscillations in the autocorrelation functions e.g. Figs. 4 and 5 
correspond to sectors [BCC) C ^ 3 (y = 2/3 ) and [BBC] C ^ 3 (y = 1/3 ), with IS's composed of repeating units BCC and 
BBC respectively. The time periods for such periodic sectors are readily calculated, since the oscillations are due to 
the alternations of l's and O's arriving at a site. 

Steady state construction. In this new scenario the initial steady configurations were prepared by means of a slight 
variation of the generating procedure described for null strings (Section 3). Consider for instance the sector [BCC] C ^ 3 . 
Then, we can identify two types of 'single site holes' B = 10 and C = 0, within an ASEP ring of L/2 + C/4 sites 
over which (L — C)/2 'monomers' (later on recast as dimers A ~ 11), are randomly adsorbed. In going backwards 
from ASEP to DDRD, we first christen the ASEP vacancies as B or C, in the order in which they occur in the 
IS. While reconstructing the DDRD configuration, we expand out the B's to be 10's and let the C's be O's. The 
resulting configuration is updated by our usual stochastic rules. This was carried out on DDRD rings of 1.2 x 10 5 
sites and averaged typically over 3 x 10 4 independent histories. Particle densities, p2 = 1 — ■§ and p\ = 1 — ^ , on 
even and odd sublattices are preserved throughout [x = (1 — C/L)/(l + C/2L)], whereas both autocorrelations and 
currents were separately computed in each sublattice. A similar SS construction and numerical considerations apply 
to the [BBC] C / 3 sector, with an evolution which now takes place in two equivalent sublattices of density p = 1 — |^ 
[x= (l-£/i)/(l + £/5L)]. 

Sublattice currents and wave velocities . For sectors such as the above, with periodic IS's, a full discussion is pos- 
sible for currents, compensating conditions, and kinematic waves with velocities given by a density derivative of the 
appropriate current in the DDRD model. The compensating condition W — U = then corresponds to the vanish- 
ing of the kinematic wave velocity given by the current derivative in the DDRD process, at which point we expect 
slower-than-exponential decays of autocorrelations. 

In these periodic sectors one can work out analytically the current on each sublattice using the DDRD- ASEP 
correspondence given in Section 2, and hence obtain the critical x c and the velocities (details can be found in [13|). 
If z\ and Z2 — 1 — z\ are the proportions of zeros on sublattices 1,2, then the current associated with the movement 
of zeros on each sublattice is 

J i = 2xz i (l-x)/[2-y(l-x)} , i=l,2. (6) 

The corresponding DDRD sublattice kinematic wave velocities Vi are obtained from Vi = d Pi Ji where pi are sublattice 
particle densities (of l's) in the DDRD. So V\,Vi vanish at a common value, x c , of the ASEP density x given by 
d x [x(l — x)/(2 — y(l — x))] = . This is a quadratic equation for x, whose root lying between and 1 is 

x c = W2(2-y)-(2-y)]/y. (7) 

Thus, x c is a monotonic decreasing function of y in the range < y < 1, varying between 0.5 and \[2— 1 . As a check 
on the result for x c , for the null sector we have y — 1, which implies x c = \[2 — 1. 

It is easy to obtain the Vi from the Ji by relating the pi to the particle density p(x) in the full DDRD (from both 
sublattices), and thence to x (in terms of which we have Ji). The result is 

Vi = 2- 4x -2/(1 -x) 2 , (8) 

independent of sublattice label i, implying that the two sublattice velocities are the same, as well as the x c at which 
they vanish. For comparison with MC results [Eq.© ], we can obtain similarly, in any sector with a periodic IS, the 
current on sublattice i of the DDRD in terms of the particle density pi on sublattice i: 

Ji = (1 - Pi ) [2z, - (2 - y) (1 - Pi)]/[2zi +y{l- Pl )\ . (9) 

The analytic results for [BCC} C ^ 3 and [BBC] C ^ 3 , obtained by inserting their respective y and Zi, are the full lines 
in the insets in Figs. [3] and [5] respectively (skewed parabolae as functions of sublattice densities). Like the null string 
currents discussed in Section 3 (inset of Fig.[T]), they agree closely with the data points resulting from the simulations. 
The sublattice autocorrelation functions shown in the corresponding main panels are at the densities of the maxima 
of the respective sublattice currents, and each show a t~ 2 / 3 power law long-time decay, verifying in turn that they are 
at the critical density. Moreover, in analogy to null string cases, it was found that on departing from these critical 
conditions the observed power law decay changes abruptly to exponential. 
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V. SPATIAL CORRELATIONS IN NULL SECTORS 



In the null sector [0] c=JVc with an even number of particles, the correlation function in the dimer problem can be 
evaluated exactly (below) using the isomorphism to non-reconstituting dimers together with the DE representation 
introduced in Section 3. Since the presence of a D particle at site i implies and is implied by the presence of an E 
particle at site (i + 1), the (unsubtracted) particle-particle correlation function can be reduced to 

(n(i) n(i + r)) = 2£(r) + £(r - 1) + £{r + 1) , (10) 

where £ (r) is the (unsubtracted) correlation function obtained from all configurations having an E particle at each of 
sites i and i + r. The weight of each such configuration can be found by mapping it to the corresponding configuration 
in the equivalent ASEP. Allowing for the weights of the £"s at i and i + r and those of m dimers and (r — 2 — 2m) holes 
between, and for multiplicities, and summing over m (from to the integer part of (r — 2)/2 ) provides the required 
£ (r). After reductions it becomes 

£(r) = [x 2 /(l + x) 2 ][l-(-xy- 1 ], (11) 

(for details, see [HI). Inserting into Eq. fTU)) above, and subtracting (n) 2 = [2x/(l + x)] 2 provides the subtracted pair 
correlation function C(r). This is oscillatory, and has a decaying envelope with correlation length £ = [ln(l/x)] _1 . 
Fig. [6] shows MC results for the (subtracted) spatial correlations C(r), in the null sector for various x. The actual 
C(r)'s show the predicted oscillation (— l) r , and the correlation lengths agree closely with the theory. 



VI. CONCLUDING DISCUSSION 



This investigation of the simplest model combining strongly broken ergodicity with driving has uncovered a wealth 
of properties, some expected and others not foreseen. In particular the model shows an interesting interplay between 
the collective (and scaling) kinetics of the basic driven model (the ASEP) and effects of the invariant IS which 
characterises the sector in the DDRD. Monte Carlo simulations shows man of these features in striking fashion. 

Despite the richness and complexity of the combined model, much has been understood analytically, for sectors 
with null or periodic IS's. The scaling behaviour seen so far appears to be in the universality class of the ASEP. The 
relationships of nonuniversal variables, like currents, densities (including critical densities), and velocities etc., have 
been largely understood. 

There are however certain important generalisations which have not been discussed here, which we now briefly 
mention: additional special sectors, alternative boundary conditions, and fc-mer generalisations. 

We can expect, following the discussion of Ref. [6( for the DRD, new effects in sectors having IS's with various 
types of structural correlation. Provided critical conditions (i£ c 's) exist in which autocorrelations show power law time 
decay, the power law for such sectors need no longer be related simply to the ASEP exponents. One would then have 
different decay patterns, according to sector ('dynamic diversity'). The combination of the methods developed here 
with those of Ref. [6[ might well allow analytic treatment of such effects. 

Changing from periodic boundary conditions to open ones with boundary injection could strongly alter the be- 
haviour. For the ASEP, that change produces non-trivial spatial correlations and dynamics and, most important of 
all, a steady state nonequilibrium phase transition. Density profiles, correlations, and dynamics differ in the different 
phases. We would expect all these features to occur in the corresponding boundary-driven version of the DDRD, in 
cases (like injection/ejection of dimers) so long as the injection and ejection processes keep the IS in tact. Again, the 
detailed properties would vary from sector to sector, and the methods used here should prove useful in this generalised 
situation. 

Finally, generalisations from reconstituting dimers to fc-mers, deserve consideration. As with dimers, the problem 
of nonreconstituting /c-mers which is of current interest lolflil ] is contained in null sectors of the general problem. The 
sector- wise mapping to the ASEP still holds, but the number of sectors is much larger for larger k, as multiplicities in 
the types of ASEP holes increase with k. The sector- wise study of the dynamics remains to be explored systematically. 
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FIG. 1: Autocorrelation functions of null strings [0] £ with different lengths L. From bottom to top they refer to C/L = 
1/5, 3/5, 1/3 and ~ y/2 — 1 (compensating condition W ~ U = 0). The latter case is consistent with a t~ 2 ^ 3 large time decay, 
whereas the other situations give rise to exponential decays. The inset displays steady currents of null sectors with particle 
densities p = 1 — C/L , closely following the analytical current p (p — l)/(p — 2) derived from Section 3, and reaching a maximum 
near vanishing wave velocities. 
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FIG. 2: Comparison of ASEP (squares) and DDRD (triangles) autocorrelation functions respectively at x — 1/2 and x — x c 
(compensating conditions in both systems). Above ~ 20 MC steps they both decay as ~ i~ 2 ^ 3 and with very similar amplitudes. 
The short time behavior is shown by the inset displaying the ratio R(t) = C DDRD (t)/Cx(t) which involves the wheeling effect 
referred to in Section 3. 
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FIG. 3: Data collapse of autocorrelations in sector [0] £ for non-critical regimes. They refer respectively to C/L = 3/5 
(squares), 1/5 (circles) and, 1/3 (triangles). The curvature of the ASEP scaling function g" (solid line, Ref. [HJ]), evidences 
that the collapse achieved through the fitting parameters a, 6 is yet very far from asymptotia. The latter is only reached above 
s 3 = at •> 200, as hinted by the actual behavior of g" exhibited in the inset. Nevertheless, the resulting values of a and b are 
understandable in terms of the analysis of g" given in Refs. flgjllif along with the wheeling velocity [Eq.([2| ]. For displaying 
purposes, all early time data (at < 1/2) have been pruned. 
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FIG. 4: Autocorrelation functions of critical BCC sector (C c /L ~ v6 — 2), for both even and odd sublattices (dashed and 
solid lines respectively). In both cases, the asymptotic behavior is consistent with a t~ 2 ^ 3 relaxation. As explained in the text, 
the early time oscillations are due to the periodicity of the four IS characters [1000]. The inset displays even and odd sublattice 
currents (upper and lower curves respectively), as functions of their corresponding densities p 2 — 1 — ~, pi = 1 — The 
data are accurately described by J2 = 2p2(l — P2)/(3 — P2) and Ji = (1 — pi)(2pi — l)/(2 — pi) , in agreement with Eq. ([9}. 
The critical condition of the main panel corresponds to maximum currents near p2 = 3 — \/6 , and pi = 2 — ^/3/2 , namely at 
vanishing wave velocities. 
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FIG. 5: Autocorrelations of critical BBC string with length density C c /L ~ \/30 — 5 . The large time behavior follows closely 
a i _2//3 power law decay (denoted by the rightmost lower dashed line). As in Fig. [4] the initial oscillations can be accounted 
for by the periodicity of the IS elements [10100]. The inset contains the steady currents (equivalent in both sublattices), which 
follow entirely their analytical counterparts J — (p — l)(5p — 2)/(p — 4) , in Eq. ((9]). The wave velocity vanishes at the current 
maximum, on approaching the main panel regime at p = 4 — 3y/ 6/5 . 
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FIG. 6: Spatial pair correlations of sector [0] £ for C/L = 1/2 (triangles), 1/3 (circles), 1/5 (squares), and 1/10 (rhomboids). 
The solid lines are fitted with slopes (inverse correlation length) = In ( ^z^j, ^ > m agreement with Eqs. (|10[l and (|lip . The 
actual oscillations of C(r) also follow the behavior predicted in Section 5. 



